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Introduction 


In  a  recent  Binary  Systems  report  [1] ,  the  difficulty  of 
measuring  the  turbulent  bounary  layer  (TBL)  wall-pressure  wavenumber 
spectrum  has  been  detailed.  It  was  shown  that  to  achieve  desired 
measurement  accuracy,  the  gain  and  phase  of  the  acoustic  elements 
must  be  very  carefully  controlled  if  conventional  measurement 
techniques  are  to  be  used.  In  fact,  the  requirements  are  so 
stringent  as  to  cast  doubt  on  the  feasibility  of  the  measurement 
program.  For  subsequent  comparison  with  the  technique  to  be 
proposed  in  this  report,  we  repeat  some  of  the  earlier  designed 
figures  here.  It  was  shown  that  the  required  array  has  40  1.2-inch 
elements  butting  against  each  other.  The  tolerance  is  0.058  dB  in 
gain,  0.3°  in  phase,  and  1.73  mil  in  phase-center  location.  Even 
with  such  tight  tolerances,  the  array  can  measure  the  wavenumber 
spectrum  only  for  values  of  (k/2ir)  between  -4.0  and  +4.0  cycle/ft 
with  3-dB  accuracy  at  UJO  Hz.  These  results  are  based  on  the 
assumption  that  the  TBL  wavenumber  spectrum  is  given  by  the  Chase 
model  [2] .  Such  a  spectrum  is  shown  in  Figure  1  where  the  frequency 
is  100  Hz  and  the  free-stream  velocity  is  25  ft/sec.  This  spectrum 
does  not  include  the  effect  of  finite  element  width.  For  thin 
elements,  the  effect  will  be  small  anyway.  Notice  the  power  levels 
at  the  convection  peak  and  at  0  cycle/ft.  The  difference  exceeds 
55  dB.  It  is  extremely  difficult  to  control  the  sidelobes  of  the 
array  factor  using  conventional  spectral  measurement  technique  so 
that  the  leakage  everywhere  is  below  the  actual  signal  level. 
Although  the  Chase  model  does  not  necessarily  represent  the  true 
TBL  spectrum,  the  measurement  technique  to  be  used  must  be  able 
to  handle  the  eventuality  that  it  is  the  correct  one.  Since  it  is 
the  more  difficult  of  the  two  models  being  suggested,  the  other 
being  the  Gardner  model  [3],  the  Chase  model  is  being  used  as  the 
design  target.  Any  technique  which  can  handle  the  Chase  model 
will  handle  the  Gardner  model  also. 

For  convenience,  the  functional  forms  of  the  two  models 
are  repeated  here.  They  are  taken  from  Appendix  A  of  Reference  [1] . 


1 


Chase  Model 


P  (k,w)  = 
o 


=t°V 


K2  Kt“5 


0.063 


v  =  0.04U 

* 


2  2 
K  =  K  + 


uj-U  k 
C 


U  =  0.68U 
c 


A  =  1.08a 


h  =  3.7 

Throughout  this  report,  k  is  wavenumber;  oi  is  angular  frequency; 

U  is  free-stream  velocity;  U  is  convection  velocity;  'a'  is 
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cylinder  radius;  p  is  fluid  density  (2  slug/ft  in  our  discussion); 
v*  is  shear  velocity.  The  spectral  density  Po(k,uj)  is  defined  so 
that  its  integral  over  both  positive  and  negative  K  and  u>  is  the 
expected  square  of  the  circumferentially  averaged  pressure. 

Gardner  Model 
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Ng  is  called  the  Strouhal  number.  Each  model  presents  the 
spectral  density  PQ(k,  w)  in  terms  of  p,  U,  and  a;  a  set  of 
fixed  numerical  parameters;  and  a  function  of  Ng  and  ka. 

In  the  discussion  to  follow  the  emphasis  is  on  the 
Chase  model.  As  seen,  the  Chase  model  is  much  more  sharply  peaked 
and  hence,  more  difficult  to  measure  for  the  reason  given  previously. 

The  spectral  estimation  technique  being  proposed  is  called 
the  maximum  likelihood  estimation  (MLE)  technique.  A  detailed 
derivation  of  this  technique  will  not  be  given  in  this  report.  The 
relevant  equations  are  given  in  Section  2.  An  intuitive  explanation 
of  the  merits  of  this  procedure  will  be  attempted.  In  Section  3, 
the  inherent  potential  of  MLE  to  measure  the  TBL  spectrum  will  be 
investigated.  Tradeoffs  of  the  number  of  elements,  element  spacing, 
and  element  size  will  be  discussed.  In  Section  4,  the  measurement 
error  induced  by  imperfect  gain  and  phase  control  will  be  considered 
It  will  be  shown  that  the  MLE  technique  alone  is  still  inadequate  for 
measuring  the  Chase  spectrum  but  a  combination  of  MLE  and  a  more 
conventional  Barlett  window  spectral  estimation  technique  seems 
quite  promising.  In  Section  5,  a  practical  implementation 

problem  of  the  MLE  will  be  considered.  The  MLE  is  based  on  the 
correlation  function  of  the  random  process  being  measured.  In 
practice  this  function  is  obtained  from  time  or  space  averaging 
which  will  not  be  perfect  but  the  effect  of  the  imperfection  on 
the  estimation  will  be  shown  to  be  small. 


In  the  concluding  section,  we  summarize  the  results 
which  have  been  obtained  so  far.  We  also  speculate  on  the 
possibility  of  making  the  element  spacing  nonuniform  to  achieve 
the  same  measurement  accuracy  with  fewer  elements. 
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2.  Maximum- Likelihood  Wavenumber  Estimation 

The  MLE  is  first  proposed  by  Capon  14]  in  1969.  A 
more  concise  discussion  of  this  technique  is  given  by  Lacoss  [5] 
in  1971.  Our  discussion  is  more  akin  to  the  latter. 


Suppose  there  is  a  receiving  array  of  N  equally  spaced 
elements  whose  outputs  are  to  be  combined  linearly  into  one 
single  sum  such  that  certain  objectives  are  fulfilled.  Let  x^ 
denote  the  outputs  of  the  elements.  The  combined  sum  z  is 
given  by 


ax 
l  l 


(3) 


In  vector  notation. 


a  x 


i  4  i 


where  both  a  and  x  are  row  vectors  of  dimension  N.  The  symbol 
’ *'  denotes  complex  conjugate  and  '+'  denotes  complex  conjugate 
transpose.  In  conventional  array  design,  the  vector  a  introduces 
a  sequence  of  delays  or  phase  rotations  such  that  the  element 
outputs  due  to  a  signal  from  a  desired  direction  will  add 
coherently.  The  undesired  'noise'  from  all  other  directions  is 
controlled  by  an  amplitude  taper  across  the  aperture  at  the 
expense  of  reduced  array  gain  and  increased  beamwidth. 
Specifically,  the  desiun  does  not  take  advantage  of  any  infcrma- 
t  ion  concerning  the  environment  even  if  it  is  .available.  On  the 
other  hand,  the  MLE  technique  measures  the  environment  and  takes 
advantage  of  it.  It  obtains  the  environmental  information  by 
formina  a  correlation  matrix  based  on  the  element  outputs.  This 
matrix  R  is  defined  by 


R 
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where  E[x]  is  the  expected  value  of  x.  It  then  derives  a  w sight 
vector  a  which  holds  the  gain  in  the  desired  direction  to  unity 
and  minimizes  the  power  from  all  other  directions.  Mathematically, 


if  S  ■  ( 1 ,  e- 


i ( N- 1 ) kA 

eJ  ) ,  where  A  is  the  element  spacing, 


represents  the  steering  vector  of  a  desired  wavenumber  k,  this 
procedure  of  finding  a  is  equivalent  to  minimizing  a '  given  by 

2  t  t 

o*Elz.z]*aRa  ( to ) 

subjected  to  the  constraint 


aSf  *  1 


The  optimization  can  be  achieved  usino  Lagrange  multipliers. 
Without  going  into  detail,  the  solution  of  a  is 


*  _  i  - 1 

S  R  SR  S 


and  the  estimated  power  at  wavenumber  k  is  aiven  by 

PL  =  (SlT1  s'1-)"1  (9) 

Intuitively,  the  MLE  would  be  superior  because  the 
environment  is  usually  inhomocreneous .  As  in  the  case  of  TBL  noise, 
the  disturbance  coming  from  certain  wavenumber  region(s)  is  often 
stronger  than  others.  A  better  receiving  array  should  obviously  be 
configured  in  such  a  manner  so  that  the  nulls  of  the  array  should 
be  placed  in  those  wavenumbers  where  the  interference  is  the 
strongest.  In  fact,  if  an  interference  is  close  to  one  side  of 
the  desired  direction,  the  peak  of  the  beam  should  be  pointed  to 
the  other  side  with  a  null  on  the  interference.  In  such  an  event, 
the  siqnal  uain  will  decrease  but  the  overall  sional-to- inter ference 
ratio  will  actually  improve. 
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In  the  formulation  above ,  the  assumption  that  the 
receiving  elements  be  equally  spaced  is  not  important.  Any 
spacing  can  be  used  as  long  as  their  exact  locations  are  known. 
As  shown  later,  it  is  rather  important  to  know  these  locations 
precisely  for  the  MLE  to  be  effective.  Any  random  error  in 
location,  gain,  or  phase  control  could  degrade  the  performance 
seriously  in  certain  wavenumber  regions.  In  order  to  overcome 
these  drawbacks,  there  will  be  a  need  to  use  the  more  conven¬ 
tional  technique  also.  As  explained  earlier,  the  weight  vector 
a  is  chosen  on  the  basis  of  the  steering  vector  alone.  In  this 
case,  the  estimated  power  is  given  by 

P_  =  -K  (SRSf)  (10) 

N 

where  S  is  again  the  steering  vector.  There  is  no  amplitude 
taper  across  the  aperture.  This  uniform  weighting  scheme  is 
sometimes  called  the  Sarlett  window.  In  (10) ,  N  is  the  number 
of  sensor  elements. 
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Inherent  Measurement  Capability  of  the  MLE 

In  this  section,  the  applicability  of  the  MLE  to  this 
measurement  program  is  explored.  Imperfections  in  the  elements 
are  ignored.  In  fact,  it  will  be  assumed  that  the  exact  spatial 
correlation  matrix  R  is  known.  Although  it  is  not  a  necessary 
condition,  it  will  be  assumed  that  the  elements  are  uniformly 
spaced.  In  particular,  we  shall  investigate  the  performance  as 
functions  of  element  spacing,  element  width,  and  the  total 
number  of  elements. 

In  all  cases  considered,  the  TBL  random  process  is 
assumed  to  follow  the  Chase  model  as  given  by  Eq.  (1) .  If  elements 
of  finite  width  are  used,  the  spectrum  will  be  modified  by  a  sin(x)/(x) 
multiplicative  factor  where  x  is  given  by 

x  =  wk/2  (11) 

and  w  is  the  width  of  the  element. 

3.1  Infinitely  Thin  Elements 

The  first  case  considered  deal  with  sensor  elements 
without  width,  i.e.,  w=0.  The  ideal  Chase  TBL  spectrum  normalized 
to  its  peak  level  P  ,  is  shown  in  Figure  1.  In  order  to  study  the 
inherent  capability  of  the  ML  method,  the  correlation  matrix  R  is 
obtained  by  inverse  Fourier  transforming  the  Chase  spectrum  as  given 

by  Eq.  (1)  directly.  Numerical  techniques  are  used.  The  span  for 

★ 

(k/2^)  is  taken  from  -24  to  +24  cycle/ft  which  should  be  wide 
enough  to  avoid  any  spectral  foldover  problem.  A  512-point  trans¬ 
form  has  been  used  in  all  the  Fourier  computations. 

Since  the  random  process  is  assumed  to  be  spatially 
stationary,  the  correlation  function  is  a  function  only  of  the 
spacing  between  the  elements,  i.e.,  R(m,n)  =  R(m-n). 


Actually  the  upper  limit  of  the  wavenumber  span  is  2tt  times  the 
inverse  element  spacing  selected.  In  most  cases  considered  in 
this  report,  the  spacing  is  0.5  inch.  Therefore  the  upper  limit 
on  (k/2T)  is  24  cycle/foot. 
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Thus,  the  correlation  matrix  needed  in  Eq .  (9)  is  given  by  (m,n 
are  indices  of  the  elements) 


where  j  is  J-T  and  A  is  the  element  spacing.  Eqs.  (9)  and  (10)  are 
used  to  compute  the  ML  and  the  conventional  estimates  respectively. 


The  first  example  is  an  array  of  25  elements  with  0.5- 
inch  spacing.  The  estimates  are  shown  in  Figure  2.  The  dash  line 
is  the  original  ideal  spectrum,  the  solid  line  is  the  ML  estimate 
and  the  dot-dash  line  is  the  conventional  estimate  with  a  Barlett 
window.  Near  the  convection  peak  both  estimates  are  reasonably 
close  to  the  true  value.  Near  zero  wavenumber,  the  conventional 
estimate  of  log  PQ(k,u))  is  30  dB  off  while  the  ML  estimate  is  12  dB 
off.  If  we  impose  a  3-dB  criterion  of  performance,  the  ML  estimate 
is  acceptable  only  for  (k/2r)  above  1  cycle/ft. 


3.2  Elements  of  Finite  Width 

In  practice,  sensor  elements  must  have  finite  width.  Its 
effect  will  be  considered  here.  Besides,  the  width  of  the  element 
will  introduce  some  spatial  smoothing  capability  which  effectively 
suppresses  the  power  of  the  higher  wavenumbers  and  improves  the 
estimates.  Most  of  the  cases  considered  below  use  0.5-inch  elements. 
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For  them,  the  first  cutoff  frequency  of  these  elements  is  at 
+24  Hz.  As  stated  previously,  the  Chase  spectrum  at  the  element 
output  will  be  modified  by  a  sin  x/x  function.  Figure  3  shows 
such  a  modified  Chase  spectrum. 

In  order  to  provide  a  comprehensive  picture  of  the  ML 
estimates,  the  number  of  array  elements  is  varied  from  5  to  25. 
Figure  4  shows  the  estimated  spectra  of  five  different  arrays. 

Notice  that  the  improvement  from  5  to  10  elements  is  considerably 
larger  than  that  from  20  to  25.  In  fact,  even  at  the  zero  wave- 
number,  the  improvement  is  only  2  dB  in  the  last  case.  By  the 
3-dB  cirterion,  the  25-element  array  is  good  for  values  of  k/2n 
to  about  0.5  cycle/ft.  It  appears  that  to  obtain  an  estimate  no 
more  than  3  dB  from  the  actual  spectrum  everywhere  the  array  size 
may  have  to  be  doubled. 

As  compared  with  the  array  with  infinitely  thin  elements, 
the  finite  element  array  is  actually  better.  This  is  due  to  the 
averaging  effect  provided  by  the  elements.  Perhaps  one  can 
improve  performance  by  further  increasing  the  element  width. 

This  is  investigated  next. 

In  Figure  5,  MLE  spectra  for  an  array  of  25  elements 
are  shown.  The  width  is  chanced  from  0.5  inch  to  0.62  and  0.74 
inches.  The  dotted  line  shown  is  the  ideal  spectrum  with  0.5-inch 
elements.  We  have  not  included  the  ideal  spectra  for  the  wider 
elements.  It  is  observed  that  there  is  some  minor  improvement 
in  the  zero-wavenumber  region  with  wider  elements,  but  the 
improvement  is  progressively  smaller.  In  fact,  the  difference 
between  0.62  in.  and  0.74  in.  elements  is  too  small  to  be 
distinguishable  in  this  figure.  Thus,  it  is  concluded  that 
element  width  is  not  a  very  critical  factor.  There  will  be 
about  1  dB  difference  between  0.5  and  0.74  in.  Subsequent  analysis 
actually  assumes  0.5-inch  elements  and  the  number  of  elements  is  set 
at  25.  A  somewhat  shorter  array  would  have  been  acceptable  also. 


venumber/ (cycle/f t) 


Figure  5.  ML  Estimates,  w=d,  25  Elements 


4.  ML  and  Conventional  Estimates  in  the  Presence  of 

Sensor  Imperfections 

The  result  in  Section  3  indicates  that  the  ML  estimation 
technique  is  acceptable  with  ideal  elements.  Although  the  spectral- 
estimate  error  at  the  zero-wavenumber  region  is  greater  than  the  desired 
limit,  the  shape  of  the  spectrum  is  reasonably  well  preserved.  If 
it  is  only  a  matter  of  choosing  the  Chase  or  the  Gardner  model,  the 
ML  estimate  will  be  sufficiently  accurate  to  effect  a  logical 

selection.  In  this  section,  the  behavior  of  the  ML  estimator  will 

be  analyzed  in  the  presence  of  element  imperfections.  Both  analytical 
and  empirical  techniques  will  be  used.  It  will  be  shown  that  the 
ML  estimate  is  quite  sensitive  to  errors  near  the  convection  peak 
but  is  quite  stable  near  the  null  region.  On  the  other  hand,  the 

conventional  Barlett  estimate  is  shown  to  be  more  stable  near  the 

convection  peak.  Together,  the  two  estimates  will  provide  adequate 
estimates  everywhere. 

Sensor  errors  are  usually  unavoidable  but  in  many 
situations  they  are  not  difficult  to  handle.  The  common  procedure 
to  overcome  these  errors  is  to  calibrate  the  sensors  by  means  of 
known  signals.  In  the  case  being  investigated,  it  is  impossible 
to  generate  a  pressure  field  of  known  frequency-wavenumber  spectrum 
to  test  the  experimental  array.  Therefore,  the  measurement 
technique  to  be  adopted  must  be  reasonably  insensitive  to  element 
imperfections.  There  are  three  kinds  of  errors  associated  with  a 
sensor  element:  gain,  phase,  and  element-center  placement.  The 
centroid  of  an  element's  sensitivity  may  not  coincide  with  the 
physical  center  of  the  element.  Any  deviation  from  the  ideal 
location  will  affect  the  estimator  performance.  If  each  deviation 
is  known,  compensation  can  of  course  be  made  and  there  will  be 
negligible  degradations.  If  the  individual  deviations  are  unknown, 
an  estimate  of  the  degree  of  degradation  in  terms  of  a  statistical 
measure  of  the  deviations  is  appropriate. 
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Observe  that  in  a  narrowband  analysis,  delay  or  dis¬ 
placement  is  equivalent  to  phase  rotation.  An  element-center 
placement  error  is  thus  equivalent  to  a  phase  rotation.  A 
signal  component  of  the  form  exp  [ j (kz-ut) ]  will  yield  an  output 
of  the  form  exp  [j  (kd-uit)  ]  from  an  element  at  d  and  exp  [  j  (kd+k6-wt)  ] 
from  an  element  at  d+5.  The  ratio  of  the  second  output  to  the 
first  is  exp  [ j k 6 ] .  Thus,  the  phase  error  equivalent  to  an 
element-center  placement  error  5  is  k5.  This  relationship  will 
be  useful  in  specifying  the  array  design.  For  the  moment  this 
relationship  allows  us  to  treat  element-center  placement  error 
as  if  it  were  just  phase  error. 

Define  an  error  matrix  T  such  that 


T  = 


ale 


-301 


*2e 


-je2 


o  .  . 


o  .  . 


Vs 


“j0N 


where  a  is  the  gain  of  the  element  and  9  is  the  phase  error 
n  n 

of  the  element.  Ideally,  an  should  all  be  unity  and  0^  should 
all  be  zero.  Let  x  be  the  input  in  the  absence  of  error.  The 
actual  received  output  will  be 

Y_  =  XT  (14) 


The  correlation  matrix  R 

y 


is  given  by 


R 


y 


il  1 


(15) 
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We  are  assuming  that  all  random  processes  have  zero  mean.  Sub¬ 
stituting  (1-1)  into  (15),  we  have 


R  -  E[(x  T)t  (x  T  )  * ) 


•Y.  T  ,  T  *  -r,  * 

=  T  E[x  x  ]  T 


-  TT  r  T  * 
x 


Substituting  (16)  into  (9),  the  maximum  likelihood  wavenumber 
spectral  estimate  is 

PL  =  <Sx^  Rx  T*)_1  Sx+)_1 


(S  (TT)-1  r  -1  T*"1  s  +)"1  (17) 


(S  R  ~1  S  + ) ~1 
y  x  y  ' 


where 


SV  -  V11’-1  *  S« 


Since  I  is  a  diagonal  matrix,  rT=r  and  the  inverse  of  T  is 
readily  obtainable.  In  fact. 


IS 


Eq.  (17)  is  important  because  it  permits  us  to  consider 
element  errors  as  errors  in  the  steering  vector  S.  In  other 
words,  it  may  be  assumed  that  the  array  is  perfect  but  the 
steering  is  wrong. 

In  the  case  of  conventional  estimates,  from  Eq.  (10), 

we  have 


Again,  we  can  treat  errors  in  the  elements  as  errors  in  the  steering 

vector  although  there  is  a  difference  between  S  and  S  ,  namely, 

- 1  i  2 

the  error  matrices  are  1  and  T  respectively. 


On  the  basis  of  the  above  discussion,  we  have  transformed 

an  imoerfect-element  Droblem  into  an  imoerfect-steer ina  Droblem. 

If  the  errors  are  small,  deviations  in  the  estimates  can  be  obtained 

analytically.  To  do  so,  we  employ  a  well  known  theorem  in 

functional  analysis  which  states  that  if  a  function  is  analytic 

near  x  ,  then 
o 
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r 


f(x  +6)  =  £  ( x  )  + 
o  o 


f'(x  )  f"(x  ) 

- -  5  +  - — 


5  + . 


(21) 


i: 


2\ 


Higher  order  terras  can  be  ignored  if  5  is  small.  In  the  present 
context,  x  is  either  gain  or  phase  error.  f(x)  is  either  P 

Li 

or  P  .  Since  it  is  assumed  that  the  errors  have  mean  zero, 


E  [  f  ( xQ+  -  f  (xQ)  ] 


f"  (x) 
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Eli*] 


(22) 


E  [5  z ]  is  the  variance  of  the  errors.  In  the  case  being  considered 

the  quantities  represented  by  f(x),  namely  P  or  P  ,  are  functions 

of  many  variables.  If  the  errors  are  basically  small,  there  will 

be  no  interaction  between  the  various  error  terms.  The  total 

degradation  is  a  linearly  weighted  sum  of  all  the  errors.  In  other 

words,  the  deviations  in  PT  or  P  are  given  by 

Lj  c 


ipL  ■  i  “s2 


E? 

*—  n 


■Pt/3  S 
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V1 


(23) 


1  2 
+  2°0 


*-  n 


J 


AP  =  5  a  2ry^32P  /  38  2_l 
c  2  8  i  L-,  c  n  j 


1  2 
+  2  a0 


E}2v3»„  j 


0  =0 

n  n 


6n=1 


0  =0 

n  n 


(24) 


where  it  has  been  assumed  the  variances  of  the  sensor  errors  are 
the  same  for  all  elements. 
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In  practice,  it  is  more  convenient  to  consider  the 
ratio  AP^/P^  or  Apc/pc'  sut>sequent  discussions,  it  is  these 

quantities  that  are  being  computed. 


The  differentiation  of  PT  and  P  is  carried  out  in 

L  C 


Appendix  A.  Defining  Q  by 


Q  =  S  R  1S  + 


b  s  s 
nm  n  m 


S 2p^  /  50^  2  and  3 ^P^  /  36 ?  2  are  given  by  ** 


—  =  2Q~3 


— 

HE 


j6£  “j0n  1  2 
Im (b.  e  1  e  ) 

x,n 


-2  ^  D0,  ~jen 

+  20  E  Se(bine  e  1 

n^  ?, 


—  =  2Q~ 3 
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'3  V 


■2q'2  b, 


The  notation  Re(x)  and  Im(x)  denote  the  real  and  imaginary 
part  of  x  respectively. 


r 


2 


are  given  by 


and  /  ^0  and  l^P^  /  36  ^ 


2  3  Pc  36  -36 

(Nz)  - =  -2  Re  ( r ,  e  e  n) 

L,  ln 
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(28) 


(N“) 


i.  i 


(29) 


where  r  are  the  elements  of  R  rather  than  R  ^  as  in  (25) . 
nm 

Substituting  the  appropriate  derivatives  into  (23)  and  (24), 
we  obtain  the  desired  analytical  expression  relating  gain  or 
phase  (placement)  errors  to  the  spectral  measurement  error. 

Note  that  the  derivatives  are  to  be  evaluated  at  the  desired  3 

n 

and  0  .  We  shall  term  these  derivatives  appearing  in  (23)  and 
(24),  the  error-sensitivity  coefficients.  If  they  are  large, 

small  errors  in  the  elements  will  cause  large  degradations.  In 

2  2  —  3  —  2 

Eg.  (26)  we  observed  that  3  P./&9.  is  dependent  on  Q  and  Q  , 

3  v. 

which  are  equal  to  P  and  P  2.  Therefore,  it  is  obvious  that 

Li  L» 

the  error-sensitivity  coefficients  will  be  large  in  the  reaion 
where  P^  is  large.  Table  la,b  is  a  summary  of  these  coefficients 
for  various  wavenumbers  assuming  the  Chase  spectrum  is  being 
measured  with  a  25-element  array  with  0.5-inch  elements  and  0.3- 
inch  spacing  between  elements.  As  seen  P  is  the  least  sensitive 

i-i 

near  zero  wavenumber  and  10,000  times  as  sensitive  near  the  con¬ 
vection  oeak.  The  sense  of  error  is  such  that  the  peak  level 
will  be  reduced  by  large  amounts.  P  ,  on  the  other  hand,  is 
fairly  insensitive  to  element  errors  in  all  regions  of  interest 
but,  relatively,  it  is  least  sensitive  near  the  convection  peak. 
Figure  6  is  the  result  of  a  sample  simulation.  The  phase  error 
is  assumed  to  be  Gaussian.  When  we  compare  this  to  Figure  6, 
it  confirms  the  conclusions  obtained  by  the  analytical  means. 


"3  •> 


The  analysis  assumes  small  errors.  An  error  less  than  3-dB  in  the  spectral  estimation  is 

deemed  acceptable. 


Note  that  the  errors  in  (23)  and  (24)  are  additive. 


From  this  error  analysis,  it  is  clear  that  relying  on 
maximum  likelihood  estimate  alone  will  impose  unrealistic  require¬ 
ments  on  the  sensors.  But  there  is  no  reason  why  one  cannot 
combine  the  two  estimates,  ML  and  conventional,  to  obtain  a 
very  respectable  estimate  with  a  much  more  practical  design. 

In  particular,  the  combined  estimate  around  the  convection  peak 
and  around  zero  wavenumber  will  be  almost  unaffected  by  sensor 
errors.  In  the  intermediate  regions,  some  interpolation  seems 

desirable . 
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5.  The  Effect  of  Imperfect  Correlation  Function 

In  the  discussion  up  to  this  point,  it  has  been  assumed 

that  an  exact  copy  of  the  correlation  matrix  R  is  available.  In 

practice,  this  matrix  must  be  obtained  from  measurements.  The 

elements  r  are  given  by 
nm 


r 

nm 


M 

E 

i  =  l 


x  ( i ) x  (i) 
n  m 


where  the  index  i  denotes  the 


.  th 
1 


time 


sample . 


(30) 


How  large  must  M  be  to  obtain  good  estimates?  An 

analytical  answer  is  believed  to  be  obtainable  but  rather  difficult. 

However,  it  can  be  shown  relatively  simply  that  M  must  be  greater 

than  N  which  is  the  number  of  elements  in  the  array.  If  this 

constraint  is  not  satisfied,  the  inversion  process  required  in 

computing  P  is  unstable  because  the  rank  of  R  would  be  less  than 
Lj 

N.  In  the  simulations  done  to  date,  M  is  usually  set  at  N+5.  Results 
obtained  are  quite  consistent  and  satisfactory.  An  example  is 
shown  in  Figure  7  for  a  25-element  array.  For  comparison,  the 
ideal  ML  estimate  is  also  shown  in  dotted  line- 
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6.  Summary 

In  conclusion,  we  have  shown  that  maximum  likelihood 
wavenumber  estimate  can  be  used  effectively  in  the  measurement 
of  flow  noise.  In  the  presence  of  sensor  error,  it  is  necessary 
to  use  a  combined  technique  to  estimate  the  power  in  different 
wavenumber  regions.  With  this  technique,  the  acceptable  phase 
error  is  on  the  order  of  3°  (including  phase-center  placement 
error,  see  Eq.  (15)) and  the  acceptable  amplitude  range  is 
between  1+0.06.  In  computing  the  correlation  matrix  for  a 
25-element  array,  an  average  of  thirty  time  samples  may  be  used. 
Briefly,  we  summarize  the  design  parameters  in  Table  2. 

Table  2.  Measurement  Array  Design 


Number  of  Elements  25 

Element  Width  0.5  inch  ! 

Element  Spacing  0.5  inch  ! 

Normalized  Amplitude  Range  1  +  0.06 

Phase  Error  1.5° 

Element-Center  Error  3  mils 


These  parameters  are  probably  more  stringent  than  necessary.  The 
analysis  on  which  these  figures  are  based  tends  to  be  pessimistic. 
As  compared  with  the  figures  given  in  Reference  [11,  the 
tolerances  are  larger  by  a  factor  of  10. 

With  this  design,  the  measurements  should  be  within 
3  dB  of  the  true  spectral  level  everywhere  where  k/27t  is  above 
0.5  cycle/ft.  Below  this  value,  the  error  is  larger  and,  at 
zero  wavenumber,  the  expected  error  is  about  6  dB. 
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In  designing  an  array,  it  is  observed  that  a  long 
array  has  better  resolution  and  better  sidelobe  suppression. 

It  is  also  observed  that  large  gaps  between  sensor  element 
could  cause  aliasing  error.  So,  if  we  are  confined  to  a 
uniform  array,  we  must  increase  the  number  of  elements  to 
obtain  better  performance  beyond  a  certain  limit.  In  Eq.  (9) , 
it  is  clear  that  the  MLE  can  be  applied  to  arrays  with  any 
type  of  element  spacing.  As  long  as  the  proper  steering  vector 
is  employed,  the  best  estimate  for  any  configuration  can  be 
obtained  just  as  easily  as  that  for  a  uniform  array.  Perhaps, 
with  nonuniform  spacing,  the  MLE  estimate  can  be  improved 
without  increasing  the  number  of  elements.  This  will  be 
investigated  in  the  future. 
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Appendix  A  The  Derivatives  of  the  Spectral  Estimates,  PT  and  P 

L  c 

From  Eqs .  (9)  and  (10),  both  PT  and  P  are  seen  to 

L  C 

have  similar  forms.  To  evaluate  the  degradation  due  to  sensor 

errors,  the  first  and  second  derivatives  are  derived  in  this  appendix. 

From  Eg.  (9),  P  is  in  the  form 
Li 


where  Q  is  the  quadratic  form 

- 1  + 
Q  =  S  R  1S ' 


b  s  s 
nm  n  m 


n ,  m 


and  from  Eq.  (10),  P  has  the  form 
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( A—  1 ) 


(A- 2) 
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P 
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=  Q' 


=  S 


R  S 


-E 

n,m 


(A- 3) 


The  summations  are  from  1  to  N.  b  are  certainly  different  from 
rnm  but  both  are  independent  of  the  variables  controlling  the 
steering  vectors  {s  }  which  is  in  the  form  of 


sn  =  8n  exp  (30n)  (A-4) 

nominally  unity  and  0^  depends  on  the  pointing  direction. 
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It  will  be  assumed  that  the  errors  are  small.  Differentiating 
PL  with  respect  to  0^,  we  obtain  from  (A-l) 


3Pt  ,  3Q 

L  -  -  (Q; 


39, 


39, 
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At  the  same  time, 
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Since  Q  and  Q'  are  essentially  the  same  function,  we  need  only  to 
2  2 

find  3  Q/3  0J,  .  We  shall  hold  constant  and  set  them  equal 
to  unity.  Substituting  (A-4)  into  (A-2)  we  have 
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we  obtain 
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where  Im(x)  is  the  imaginary  part  of  x.  Differentiation  of 
time  with  respect  to  9 ^  to  obtain  3 2Q/3 9  2  yields 
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where  Re(x)  is  the  real  part  of  x.  Combining  (A-6),  (A-12)  and 
(A-13) ,  we  obtain  the  desired  result  for  32Pl/  30e  2.  Naturally, 
it  we  replace  b^m  in  (A-13)  by  r^m  ,  we  obtain  32P  / 36  2 . 

O  X, 

The  derivative  of  and  P^  with  respect  to  B^can  be 
found  in  a  similar  way.  Eq.  (A-6)  and  (A-8)  will  hold  if  3Q 

is  replaced  by  30/36^.  In  this  case  0n  are  held  constant .  e For 

convenience,  in  (A-14)  to  (A-16)  they  are  treated  as  part  of 
b  or  r 

nm  nm.  Hence, 
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The  first  derivative  of  Q  is 


'  2  \t  V  E  <b«»  6m  +  b„,<  V 


=  2ERe<btm  6m) 


(A-15) 


and  the  second  derivative  is  simply 
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Combining  (A-6),  (A-15)  and  (A-16),  the  desired  3P  2/ 36 .  is 
obtained,  i.e. 
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